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Abstract. - Quantum decay in an ac driven biased periodic potential modeling cold atoms 
in optical lattices is studied for a symmetry broken driving. For the case of fully chaotic 
classical dynamics the classical exponential decay is quantum mechanically suppressed for a 
driving frequency uj in resonance with the Bloch frequency u>b, qu = tub with integers q 
and r. Asymptotically an algebraic decay ~ i -7 is observed. For r = 1 the exponent 7 
agrees with q as predicted by non-Hermitian random matrix theory for q decay channels. The 
time dependence of the survival probability can be well described by random matrix theory. 
The frequency dependence of the survival probability shows pronounced resonance peaks with 
sub-Fourier character. 



Introduction. - Chaotic classical dynamics has quantum signatures in statistical prop- 
erties of eigenvalues, wave amplitudes and state projections [1] which can be described by 
random matrix theory. This has been demonstrated in many cases for bound systems. For 
open systems, however, such studies are rare. Only recently, the statistics of the lifetimes 
of the metastable resonance states (or the complex part T of the resonance energies) have 
been shown to be in agreement with the appropriate random matrix results. This has been 
achieved for only few physical systems, e.g. open billiards [2-4], scattering on graphs [5] and 
strongly ac driven periodic lattices [6,7] which will also be studied in the present letter. 

Open classical systems with fully chaotic intrinsic dynamics show generically an exponen- 
tial decay of the survival probability P(t) = e~ vt in the long time limit, where the decay 
rate v can be related to the Lyapunov exponent and the fractal dimension of the chaotic 
repeller [10]. Quantum mechanically, the situation is more complicated because the number 
of open decay channels depends on the Planck constant h and typically increases in the limit 
of small H. 
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The model. - The one-dimensional Hamiltonian 



H = p 2 /2 + V(x) + F(t)x , V(x + 2vr) = V(x) , 



(1) 



with space period d = 2n and ac-dc driving 



F(t) = F + F w (t) , F w (i + TJ) = F„ (t) 



(2) 



(where the time-average of F w {t) can be chosen to be zero) is known as the ac-dc Wannier- 
Stark system. Such systems have been studied by many authors in different context (see, 
e.g. the review [7]). In recent years, one could observe an increasing interest in view of the 
exploding studies on the dynamics of cold atomic gases or Bose-Einstein condensates in optical 
lattices, where the static field Fq is generated by chirping the laser frequencies or simply by 
the gravitational force. 

Note that in the scaled units used here, h depends on the system parameters and can be 
expressed as h — 4y Er/Vq where Vb is the potential depth and Er = h 2 k 2 / (2m) is the recoil 
energy [7] . In present experiments a ratio of Vb /Er ~ 500 is routinely available which yields 
h « 0.2 not far from the value K = 0.1 used in the numerical calculations below. 

It is convenient to rewrite the dynamics in a Kramers-Henneberger form [7] 



with K u (t) = f t dt' G u (t') , G u (t) = ^dt' F u {t') . We will study the potential V{x) = cos a; 
with de-symmetrized driving 



and we choose the time t as a solution of 2 sin wig = cos 2tot with the consequence that 



(e = A^/lu 2 ) is periodic. 

Chaotic dynamics. - Let us first consider the classical dynamics. A stroboscopic phase 
space plot for the non-biased case Fq = 0, e.g. for parameters u> = 1, e = 3, reveals a chaotic 
strip with few embedded regular islands. When a weak static force Fo is added, this chaotic 
strip survives but the invariant curves confining the chaotic strip are destroyed. A trajectory 
starting inside the strip will show a diffusive motion until it reaches the (former) boundary pi 
where it escapes to infinity. As a measure of the local spreading in phase space we calculate 
the stability (or monodromy) matrix M(t) ( x ) for a trajectory started at a phase space point 
(x,p). Figuren]shows the norm ||M(t)|| at time t = 6T W as a function of (x,p) for ui = 1, e = 3 
and an additional static force Fo = 0.016, a value used also in the quantum case discussed 
below for resonant driving. We clearly observe the chaotic strip where the lower boundary 
at pi Ri —4 is relatively sharp whereas the upper boundary at p2 ~ +6 is more diffusive. 
This difference is due to the fact that particles with large enough negative momentum escape 
and particles with positive momentum are reflected by the increasing linear potential and 
re-injected into the chaotic strip. A similar behavior is found for other values of Fo used in 
the quantum calculations discussed in the following. 

( 1 )For a bound system, the long time average of the norm of the stability matrix gives the Lyapunov exponent. 
Here the majority of the trajectories escape to infinity and this limit is not defined. For more information on 
phase space derealization see [8]. 



#kh - (P - F a t) 2 /2 + V(x - K u (t) ) 



(3) 



Fuj (i) = A u ( cos cot + sin 2ujt ) 



(4) 



Kcj (t) — — — (4 cos tot + sin 2uit — 4 cos ujto — sin 2wto) 



(•5) 
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Fig. 1 - Norm of the monodromy matrix as a function of the initial phase space point (x,p) for a time 
propagation over six periods (left) and time evolution of a wavefunction \if)(p,t) \ initially localized in 
the chaotic strip (right). The parameters are u) = 1, e = 3 with a static field is Fo = 0.016. 



The quantum dynamics is qualitatively different because of the Bloch oscillation [9], an 
intrinsic periodic motion with the Bloch period in the non-driven case F u (t) = 0. By con- 
structing a generalized scattering-matrix for such systems, resonance states can be defined 
and also easily computed [7], provided that the driving and the Bloch frequency lub = dFo/h 
are in resonance, 

quartos, (6) 

where r and q are coprime integers. The integer q is equal to the dimension of the scattering 
matrix, i.e. the number of open channels [7]. This offers the unique possibility to tune the 
number of decay channels by varying the external system parameters. 

Here we will explore the quantum decay dynamics for r — 1 as a function of the num- 
ber q of decay channels, starting from the case of resonant driving, to = lub, i.e. q = 1. 
Here and in the following we use h — 0.1. The quantum wavefunction \ip(p, t)\ is plotted in 
fig. ^for an initial minimum uncertainty Gaussian wavepacket with momentum width 0.28 
centered at p = in the classically chaotic region. Also here the chaotic strip can be easily 
identified: The wavefunction rapidly spreads over the strip and decays in the direction of 
negative momentum. Inside the strip we find an irregular oscillation and outside, for p < pi, 
a decay to — oo with linearly growing momentum p ~ —Fat. If the assumption of an irregular 
wavefunction inside the chaotic strip is correct, the statistics of the normalized probabilities 
s n = \ip(pn, t)\ 2 1 Yin ^(Pm*)! 2 a ^ discrete values of the momentum should follow the ran- 
dom vector model [1, 12, 13] describing in general the fluctuations of a sum of v independent 
Gaussian distributed variables (y = 1,2,4) for the orthogonal (GOE), unitary (GUE) and 
symplectic (GSE) ensemble, respectively. This model predicts approximately a \t distribu- 
tion 

Figure [5] shows a comparison of the numerical data of fig. ^ in the chaotic strip for t = 
100 T u with the distributions JJJ for v = 1,2,4. As expected, we observe an agreement 
with the GUE statistics. Moreover, this behavior is found to be independent of the ratio 
lu/ojb — 1/q. The decay, however, is very sensitive with respect to this ratio. 
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Fig. 2 - Statistical distribution of s = \ip(p, 100TL,)| 2 within the chaotic strip for q—1 (+) and q — 2 
(o) in comparison with the random matrix predictions J7J for the GOE (— ■ — ■), GUE ( — ) and the 
GSE ( ) case. Shown is the distribution VK(log a;) with x = s/s in order to increase the differences . 



Dynamics of decay. - In the following, we will explore the validity of the conjecture that 
the decay of these Wannier-Stark systems can be described by non-Hermitian random matrix 
theory for systems with q decay channels. 

According to random matrix theory, the distributions of the resonance widths V of the 
Wannier-Stark resonance states (the eigenstates of the Floquet operator) are given by [14] 



(—'\') q di 

nrr ) - y ' f"- 1 
{ s) (g-l)! * dT! 



1 - exp(-2r s ) 



2r, 



(8) 



for the circular unitary ensemble (CUE). Here q is the number of channels and T s = ttT/A 
is the width scaled by the mean level spacing A. Equation JSJ) can also be aplied to the case 
of a harmonic driving as shown in [7] if one restricts the Hamiltonian to a fixed value of the 
quasimomentum k which is a good quantum number. With the exception of the center and 
the boundaries of the Brillouin zone the time reversal symmetry is then broken. 

A random matrix formula for the time dependence of the decay of an open quantum system 
has been derived by Savin and Sokolov [15] using supersymmetry techniques. For CUE the 
case of perfect coupling is realized [7] (i.e. one has T = 1 in cq. (11) of ref. [15]) and the decay 
probability depends only on the number of channels q and the ratio r = t/Tjj where Th is 
the Heisenberg time Th = 2ttH/A: 



P (q) (t) = \ 



du 
i j\ 



dv o (Zt+u — v) 



1 1 



The integral can be expressed in terms of hypergeometric functions as 

P (9) (t) = (l+r)-«/<(r) forrgl, 



(9) 



(10) 



fc=0 



1 + T 
fc+1 



2^1 (q,k+l;k+2; 



2t f t 
■ 2 F 1 [q,k+2;k+3- 



k + 2 



1+T 



1+r 



(11) 
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/>(r) = ^ T f 2 F 1 (« ) l;, + 2; IT -) 

-(^^) 2Fl (^ 2 ^ +3; TT7) 

Asymptotically the decay is algebraic, 

P (q) (t) -> T- q /(q + 1) forr^oo, (13) 

where the exponent is equal to the number of decay channels. 

Quantum stabilization, i.e. algebraic decay was already observed in previous studies for a 
purely harmonic driving (see [16] and references therein) but with markedly different decay 
exponents. For the symmetry broken driving (@J) (and therefore with proper CUE symmetry 
contrary to the harmonic case) also good quantitative agreement with the channel number q 
is found as we will show in the following. 

The system parameters are the same as in fig. the static force Fq is tuned to give a 
rational ratio lu/lob = r/q with r = 1. A minimum uncertainty wavepacket with momentum 
width 0.28 centered at p = and x = it is propagated in time and the survival probability in 
the chaotic strip p\ < p < p2 is calculated from the momentum distribution, 

P(t)= / dp \i>(p,t)\ 2 , (14) 
J pi 

where p\ = — 5 and P2 = 7 are chosen based on fig-^ Numerically, these survival probabilities 
are calculated up to a time t — 200 T u for q = 1, 2, 3 and 4 . 

A fit of the long time behavior of the numerical data to the asymptotics Ijl3(l confirms the 
algebraic decay predicted by random matrix theory: Numerically the exponents are calculated 
as 0.999, 2.006, 3.094 and 4.087 for q = 1, 2, 3 and 4 in precise agreement with these channel 
numbers (compare the double logarithmic plot in fig.[2I)- The corresponding Heisenberg times 
are found to be Th/T u = 53.2, 36.5, 30.0 and 25.3 and decrease approximately as ~ q~ x / 2 . For 
resonant driving © with r ^ 1 the decay is also found to be asymptotically algebraic, however 



(12) 




Fig. 3 - Decay of the quantum survival probability P(t) for rational frequency ratios u/lub = 1/g 
with q = 1, 2, 3, 4 (open circles, from top to bottom) in comparison with the results given by non- 
Hermitian random matrix theory for q decay channels (solid lines). Right: double logarithmic plots 
with asymptotic linear fits. 
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Fig. 4 - Quantum survival probability Pit) for irrational frequency ratios uj/lub = I/a/2 (full curve) 
in comparison with the results from a corresponding classical ensemble (open circles). 

with exponents differing from q which is presumably due to the fact that the quasienergy 
spectrum is r-fold degenerate which introduces additional symmetries into the system. 

Figure [31 shows a logarithmic plot of the numerical results for the survival probabilities 
(I14fl as a function of t/T u up to t — 20021, for rational ratios lo/lob = l/<7 for q = 1, 2, 3 
and 4. Also shown in fig. [3] is the time dependence of the random matrix model IjlUfl using 
the values of the Heisenberg time extracted by the asymptotic fit. In order to account for the 
initial spreading of the wavepacket to equilibrate over the chaotic strip, the numerical data 
are shifted by 8 T u in each case. Good agreement is observed. The small deviations between 
the numerical data and the random matrix predictions (in particular for q — 3) are due to the 
influence of small stability islands in the chaotic strip. 

For an irrational ratio lo/lob the quantum survival probability decays exponentially, Pit) = 
e~ vt , as the classical distribution and even the exponent v agrees with the classical one. A 
comparison is given in fig. 01 As shown by Puhlmann et al. [17], the random matrix predictions 
for the survival probability can also be analyzed using a semiclassical approach. A semiclassical 
treatment of the driven Wannier-Stark system along these lines, however, deserves further 
studies. 

Resonance profiles. - It is obvious that the quantum survival probabilities P(t) depend 
very sensitively on the frequency ratio lo/lob- Asymptotically, the decay is algebraic for ratio- 
nal and exponential for irrational ratios and therefore a fractal-like spike structure develops as 
already discussed in ref. [11]. As a typical example, fig.0shows Pit) as a function of lo/lob in 
the vicinity of the resonance lo = lob for t — 300XL,. We observe a pronounced resonance peak 
which sharpens with increasing time. The maximum decays as ~ i" 1 and the almost constant 
background decays exponentially. The peak profiles can be well reproduced by the functional 
form Pfit(t) = a j y ' (lo — lob) 2 + 3(Alo) 2 + b also shown in the figure. The parameters a, b and 
Alo are functions of time. The time decay of the width Alo of the distribution was found to 
be algebraic, Aw ~ i -7 , with 7 sa 1.45. Interestingly, this is much faster than the Fourier be- 
havior where two frequencies can only be distinguished if their distance is at least of the order 
the excitation time, i.e. Alo ~ t . Similar observations of such 'sub-Fourier' resonances have 
been reported recently in an experimental study of bichromatically kicked cold atoms in an 
optical lattice [18,19]. This system also shows quantum chaotic behavior and the mechanism 
of the sub- Fourier sharpening of these resonances could be explained theoretically [20] based 
on the phenomenon of dynamical localization. 
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Fig. 5 - Survival probability P(t) as a function of lo/lub in the vicinity of the resonance uj = ujb for 
t = 3007L. 
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